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Abstract 

Fundamental aspects of inverse kinetic theories for incompressible Navier-Stokes equations con- 
cern the possibility of defining uniquely the kinetic equation underlying such models and further- 
more, the construction of a kinetic theory implying also the energy equation. The latter condition 
is consistent with the requirement that fluid fields result classical solutions of the fluid equations. 
These issues appear of potential relevance both from the mathematical viewpoint and for the phys- 
ical interpretation of the theory. In this paper we intend to prove that the non-uniqueness feature 
can be resolved by imposing suitable assumptions. These include, in particular, the requirement 
that the kinetic equation be equivalent, in a suitable sense, to a Fokker-Planck kinetic equation. 
Its Fokker-Planck coefficients are proven to be uniquely determined by means of appropriate pre- 
scriptions. In addition, as a further result, it is proven that the inverse kinetic equation satisfies 
both an entropy principle and the energy equation for the fluid fields. 
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I. INTRODUCTION 



An aspect of fluid dynamics is represented by the class of so-called inverse problems, in- 
volving the search of model kinetic theories able to yield a prescribed set of fluid equations, 
with particular reference to the continuity and Navier-Stokes (N-S) equations for Newtonian 
fluids, by means of suitable velocity-moments of an appropriate kinetic distribution func- 
tion /(r, v,t). Among such model theories, special relevance pertains to those describing, 
self-consistently, isothermal incompressible fluids described by the so-called incompressible 
Navier-Stokes equations (INSE) in the sense that the fluid equations are satisfled for arbi- 
trary fluid fields and for arbitrary initial conditions for the kinetic distribution. In particular 
the fluid flelds are in this case identifled with {p = p^, V,p} , being Vthe fluid velocity and, 
pand ^respectively, the mass density and fluid pressure, both non-negative in the closure [fl) 
of the fluid domain Q C M^, while pis always constant (condition of incompressibility) . A 
desirable feature of the theory is, however, the requirement that the relevant inverse kinetic 
equation are uniquely deflned (in particular, by imposing that a particular solution is pro- 
vided by local Maxwellian equilibria) as well as the possibility of requiring, besides INSE, also 
additional fluid equations, to be satisfled by means of suitable moment equations {extended 
INSE). An example is provided by the energy equation, i.e., the fluid equation obtained by 
taking the scalar product of the Navier-Stokes equation by the fluid velocity V 

+ V- V— +— V- Vp + — V-f-z/V-V2V = 0, (1) 

dt 2 2 po Po 

where f is the volume force density acting on the fluid element. In fact, it is well known that 
the energy equation is not satisfled by weak solutions of INSE and, as a consequence, also 
by certain numerical solutions, such as those based on weak solutions such as possibly so- 
called flnite-volume schemes. Therefore, imposing its validity for the inverse kinetic equation 
yields is a necessary condition for the validity of classical solutions for INSE. In a previous 
work , an explicit solution to INSE has been discovered based on a continuous 

inverse kinetic theory, adopting a " Vlasov " differential kinetic equation deflned by a suitable 
streaming operator L. Basic feature of this kinetic equation is that, besides yielding INSE 
as moment equations, it allows, as particular solutions, local kinetic equilibria for arbitrary 
(but suitably smooth) fluid fluids {po, V,p} . However, as pointed out in Refs. , the 
inverse kinetic equation defined in this way results parameter-dependent and hence non- 
unique, even in the case of local Maxwellian kinetic equilibria. This non-uniqueness feature 



may result as a potentially undesirable feature of the mathematical model, since it prevents 
the possible physical interpretation of the theory (in particular, of the mean-field force F) 
and may result inconvenient from the mathematical viewpoint since the free parameter may 
be chosen, for example, arbitrarily large in magnitude. It is therefore highly desirable to to 
eliminate it from the theory [4,]. 

The purpose of this paper is to present a reformulation of the problem which permits to 
cast the inverse kinetic equation in a form which results unique, thus eliminating possible 
parameter-dependences in the relevant streaming operator. Actually, the prescription of 
uniqueness on the kinetic equation is to be intended in the a suitably meaningful, i.e., to 
hold under the requirement that the relevant set of fluid equations are fulfilled identically 
by the fluid fields in the extended domain Qx I. This means that arbitrary contributions in 
the kinetic equation, which vanish identically under a such an hypothesis, can be included 
in the same kinetic equation. Consistent with the previous regularity assumption, here we 
intend to consider, in particular, the requirement that the inverse kinetic equation yields 
also the energy equation ([T]). Note that this new formulation of the inverse kinetic theory is 
also relevant for comparisons both with previous literature dealing with the determination of 
the probability distribution function (PDF) for incompressible fluids [5] and with emerging 
theoretical approaches for the determination of the PDF for small scale turbulence to model 
homogeneous and isotropic turbulence in the inertial range. As further development of the 
theory, it is shown that the streaming operator can be suitably and uniquely modified in such 
a way that the inverse kinetic equation yields the extended INSE equations, i.e., besides the 
incompressible Navier-Stokes equations also the energy equation. In particular we intend to 
prove that the mean-field force F can be uniquely defined in such a way that both kinetic 
equilibrium and moment equations yield uniquely such equations. 

II. NON-UNIQUENESS OF THE STREAMING OPERATOR 

We start recalling that the inverse kinetic equation which is assumed of the form 



[3(]. In particular, the streaming operator L is assumed to be realized by a differential 
operator of the form L{F) = §1 + v-^ + ^ ■ {F} . The vector field F (mean-field force) can 




(2) 
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be assumed of the form F = Fq + Fi, where Fq and Fi, requiring in particular that they 
depend on the minimal number of velocity moments ^, can be defined as follows 

Fo(x,t; f) = - [V-n - Vpi - f] + u • VV+z/V^V, (3) 

Po 

Fi(x,t; /) = iu Inpi+lv ■ Q-^Vp-q| + ^Vp - || . (4) 
2 [Dt pi Pi J 2pi {vf^ 2 J 

Here Q and n are respectively the relative kinetic energy flux and the pressure tensor 
Q = Jd^vuff, n = Jd^vuuf. As a consequence, both Fq and Fiare functionally de- 
pendent on the kinetic distribution function /(x,t). Supplemented with suitable initial and 
boundary conditions and subject to suitable smoothness assumptions for the kinetic distri- 
bution function /(x,t), several important consequences follow 

• the fluid fields {po,V,p} can be identified in the whole fluid domain Q with suit- 
able velocity moments (which are assumed to exist) of the kinetic distribution func- 
tion /(x,t) [or equivalent f{'x,t)], of the form MG{r,t) = f d^vG{yi,t)f{'x,t), where 
G(x, t) = 1, \,E = ^u"^, vE, uu, and u = v — V(r,t) is the relative velocity. Thus, we 
require respectively Po = J d^vf{yi,t), V(r,t) = ^/ '^'^^v/(x,t), p{r,t) = pi(r,t) — Po 
Pi{r,t) being the scalar kinetic pressure, i.e., pi{r,t) = J (ivy-/(x,t). Requiring, 
Vp(r,t) = Vpi(r,t) and pi{r,t) strictly positive, it follows that Po is an arbitrary 
strictly positive function of time, to be defined so that the physical realizability con- 
dition p{r,t) > is satisfied everywhere in x J (/ C M being generally a finite time 
interval); 

• {po, V,p} are advanced in time by means of the inverse kinetic equation Eq.([2]); 

• By appropriate choice of the mean-field force F, the moment equations can be proven 



to satisfy identically INSE, and in particular the Poisson equation 
as well the appropriate initial and boundary conditions (see Ref. 



br the fluid pressure, 

3); 



The mean-field force F results by construction function only of the velocity moments 
|po, V,pi, Q,n} , to be denoted as extended fluid fields. 

In particular, L(F) can be defined in such a way to allow that the inverse ki- 
netic equation admits, as a particular solution, the local Maxwellian distribution 
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fuij^it] V,pi) = 3 exp {— X^} .Here, the notation is standard |3|, thus X"^ = 
v^f^ = 2pi/po,pi being the kinetic pressure. 

Let us now prove that the inverse kinetic equation defined above ([2]) is non-unique, even 
in the particular case of local Maxwellian kinetic equilibria, due to the non-uniqueness in 
the definition of the mean-field force F and the streaming operator L(F). In fact, let us 
introduce the parameter-dependent vector field F{a) 

F(a) = F + au- VV-«VV-u = Fo(«) + Fi (5) 

where F = F(a = 0), a G M is arbitrary and we have denoted 



Fo(«) = Fo - aAFo = Fq^ + AiFo(a), 

AFo = u ■ VV-VV ■ u, (6) 

AiFo(a) = (1 + a)u ■ VV - aVV ■ u, 

where Fq and Fi given by Eqs.([3]),(llD. Furthermore, here we have introduced also the 
quantity AiFo(a;) to denote the parameter-dependent part of Fo(q;). In fact, it is immediate 
to prove the following elementary results: 

a) for arbitrary a; G M, the local Maxwellian distribution Jm is a particular solution of 
the inverse kinetic equation ([2]) if and only if the incompressible N-S equations are satisfied; 

b) for arbitrary a in R, the moment equations stemming from the kinetic equation 
coincide with the incompressible N-S equations; 

c) the parameter a results manifestly functionally independent of the kinetic distribution 
function /(x, t). 

The obvious consequence is that the functional form of the vector field Fq, and conse- 
quently F, which characterizes the inverse kinetic equation (j2]) is not unique. The non- 
uniqueness in the contribution Fo(a) is carried by the term aAFo which does not vanish 
even if the fiuid fields are required to satisfy identically INSE in the set Q x I. We intend 
to show in the sequel that the value of the parameter a can actually be uniquely defined by 
a suitable prescription on the streaming operator and the related mean-field force. 
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III. A UNIQUE REPRESENTATION 



To resolve the non-uniqueness feature of the functional form of the streaming operator L, 
due to this parameter dependence, let us now consider again the inverse kinetic equation ([2]). 
We intend to prove that the mean-field force F, and in particular the vector field Fo(a), can 
be given an unique representation in terms of a suitable set of fluid fields {po, V,pi, Q,n} 
defined above by introducing a symmetrization condition on the mean field force Fo(a).To 
reach this conclusion it is actually sufficient to impose that the kinetic energy flux equation 
results parameter-independent and suitably defined. Thus, let us consider the moment 
equation which corresponds the kinetic energy flux density (ji'(x, t) = v^. Requiring that 
f{^,t) is an arbitrary particular solution of the inverse kinetic equation (not necessarily 
Maxwellian) for which the corresponding moment q = / d^vv^f (kinetic energy flux vector) 
does not vanish identically, the related moment equation takes the form 

^ J civG(x,t)/ + V- J dvvG(x,t)/- 

-Jdv [Foa + AiFo(a) + Fi] ■ ^^^/- (7) 
dvf 



^G(x,t)+vVG(x,t) 



0. 



Introducing the velocity moments P2 = J dv^f, P = / dvuu^f and T = J dvuuu/, 
the kinetic energy flux equation contains contributions which depend linearly on the unde- 
termined parameter a. The contribution to the rate-of-change of q produced by AiFo(a), 
which results proportional both to the velocity gradient VV and the relative kinetic energy 
flux Q, reads 

M„(/) = -(l + a)Q-VV + aVV- Q. (8) 

In order to eliminate the indeterminacy of a, since a cannot depend on the kinetic distri- 
bution function /, a possible choice is provided by the assumption that Ma(/) takes the 
symmetrized form 

M.(/) = -^VV-Q+1q-VV, (9) 

which manifestly implies a = 1/2. Notice that the symmetrization condition can also be 
viewed as a constitutive equation for the rate-of-change of the kinetic energy flux vector. 
In this sense, it is analogous to similar symmetrized constitutive equations adopted in cus- 

n 

tomary approaches to extended thermodynamics ^]. On the other hand, Eq.([9]) implies 



= |Qx^, ^ = V X V being the vorticity field. Thus, Mq,(/) can also be interpreted 
as the rate-of-change of the kinetic energy flux vector Q produced by vorticity field C,- From 
Eq.Q it follows that Fo(a) reads 

Fo(a = h = — [V-n - Vpi - f] + ^ (u ■ VV+VV • u) + z/V'V. (10) 

Hence, the functional form of the streaming operator results uniquely determined. As a 
result of the previous considerations, it is possible to establish the following uniqueness 
theorem: 

THEOREM 1 — Uniqueness of the Vlasov streaming operator L(F) 

Let us assume that: 

1) the fluid fields {p, V,p} and volume force density f(r, V, t) belong respectively 
to the functional settings {V(r,t),p(r,t) G C^^^iTl x /), V(r,t),p(r,t) G x /)} and 
{f(r, v,t) G CW(n X /),f(r,t) G C^^'^\n x /)} ; 

2) the operator L(F), defining the inverse kinetic equation 1^, has the form of the Vlasov 
streaming operator L; 

3) the solution, /(x, t), of the inverse kinetic equation ^ exists, results suitably smooth 
in T X I and its velocity moments {po, V,pi, Q, n} define the fluid fields {po, V,p} which 
are classical solutions of INSE, together with Dirichlet boundary conditions and initial con- 
ditions. In addition, the inverse kinetic equation admits, as particular solution, the local 
Maxwellian distribution fM', 

4) the mean-field force F{a) is a function only of the extended fluid fields 
|po, V,pi, Q,n} , while the parameter a does not depend functionally on /(x, t); 

5) the vector field AiFo(a) satisfies the the symmetry condition 

Then it follows that the mean-field force F in the inverse kinetic equation ^ is uniquely 
defined in terms o/ F = Fq + Fi,where the vector fields Fq and Fi are given by Eqs. ^^) 
and 0). 

PROOF 

Let us consider first the case in which the distribution function f{^,t) coincides with 
the local Maxwellian distribution /m. In this case by definition the moments Q,n vanish 
identically while, by construction the mean mean-field force is given by F(a) [see Eq.(|5])], 
a G M being an arbitrary parameter. Let us now assume that /(x, t) is non-Maxwellian and 
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that its moment Mq,(/) defined by Eq.Q is non- vanishing. In this case the uniqueness of F 
follows from assumptions 4 and 5. In particular the parameter a is uniquely determined by 
the symmetry condition ^ in the moment Mq(/). Since by assumption a is independent of 
/(x, t) the result applies to arbitrary distribution functions (including the Maxwellian case). 
Let us now introduce the vector field F'= F+ AF, where the vector field AF is assumed to 
depend functionally on /(x, t) and defined in such a way that: 

A) the kinetic equation L(F')/(x, t) = yields an inverse kinetic theory for INSE, satis- 
fying hypotheses 1-5 of the present theorem, and in particular it produces the same moment 
equation of the inverse kinetic equation ([2]) for Giy^fy = 1, y,E = ^u^; 

B) there results identically AF(/m) = 0, i.e., AF vanishes identically in the case of a 
local Maxwellian distribution Jm- 

Let us prove that necessarily AF(/) = also for arbitrary non-Maxwellian distributions 
/ which are solutions of the inverse kinetic equation. First we notice that from A and B, due 
to hypotheses 3 and 4, it follows that AF must depend linearly on Q,n— pil- On the other 
hand, again due to assumption A the vector field AF must give a vanishing contribution 
to the moments the kinetic equation evaluated with respect to G(x,t) = l,v,E = 
Hence, in order that also F' depends only on the moments |po, V,pi, Q,n} (hypothesis 4) 
necessarily it must result AF(/) = also for arbitrary non-Maxwellian distributions /. 

IV. FULFILLMENT OF THE ENERGY EQUATION 

As a further development, let us now impose the additional requirement that the inverse 
kinetic theory yields explicitly also the energy equation ([1]). We intend to show that the 
kinetic equation fulfilling such a condition can be obtained by a unique modification of the 
mean-field force F = Fo(x,t) + Fi(x,t), in particular introducing a suitable new definition 
of the vector field Fi(x,t). The appropriate new representation is found to be 



Fi(x,t;/) 





(11) 



As a consequence, the following result holds: 
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THEOREM 2 - Inverse kinetic theory for extended INSE 

Let us require that: 

1) assumptions 1-3 of Thm.l are valid; 

2) the mean-field F is defined in terms of Fq and Fi given by Eqs. ^ and / flT]) . 
Then it follows that: 

A ) {p, V,p} are classical solutions of extended INSE in Qxl if and only if the Maxwellian 
distribution function fM is a particular solution of the inverse kinetic equation 

B) provided that the solution /(x,t) of the inverse kinetic equation ^ exists in T x I 
and results suitably summable in the velocity space U, so that the moment equations of ^ 
corresponding to the weight- functions G(x,t) = 1, v,^' = ^u^ exist, they coincide necessarily 
with extended INSE. 

C) the two representations 0) and Dl]) for Fi coincide identically 
PROOF: 

A) The proof is straightforward. In fact, recalling Thm.l, in j^, we notice that Eqs. ( fTTjl 
and (jll) manifestly coincide if and only if the energy equation ([T]) is satisfied identically, i.e., 
if the fiuid fields are solutions of extended INSE. 

B) The first two moment equations corresponding to G'(x,t) = 1, v are manifestly inde- 



lian and non-Maxwellian distributions, 
sl applies, i.e., the moment equations 



pendent of the form of Fi, both in the case of Maxwe 
i.e., (11 II) and (jlj). Hence, in such a case Thm.3 of 
yield INSE. Let us consider, in particular, the third moment equation corresponding to 
G(x,t) = 1^2, 

^Pi + V ■ Q + V ■ [Ypi] - ^ y rfvF(x,t)u/ + ^VV : n = 0. (12) 

Invoking Eqs. flTUl) and flTTl) for Fq and Fi, the previous equation reduces to piV ■ V = if 
and only if the energy equation ([T]) is satisfied. Since by construction pi > 0, this yields the 
isochoricity condition V • V = 0. 

C) Finally, since thanks to A) {p, V,p} are necessarily classical solutions of INSE, it 
follows that they fulfill necessarily also the energy equation ([1]). Hence, (jl]) and (ITT!) coincide 
identically in Fx/. 

We conclude that (JTOl) and ( JTTl) provide a new form of the inverse kinetic equation ap- 



plying also to non-Maxwellian equilibria, which results alternative to that given earlier in 
sl. The new form applies necessarily to classical solutions. Since weak solutions (and hence 
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possibly also numerical solutions) of INSE may not satisfy exactly the energy equation, the 
present inverse kinetic theory based on the new definition given above [see Eq.(jl])] for the vec- 
tor field F(x,t) provides a necessary condition for the existence of strong solutions of INSE. 
The result seems potentially relevant both from the conceptual viewpoint in mathematical 
research and for numerical applications. 



V. CONCLUSIONS 



In this paper the non- uniqueness of the definition of the inverse kinetic equation defined 
by Ellero and Tessarotto (see [sl) has been investigated, proving that the mean-field force 
F characterizing such an equation depends on an arbitrary real parameter a. To resolve 
the indeterminacy, a suitably symmetrization condition has been introduced for the kinetic 
energy flux moment equation. As a consequence, the functional form the mean-field force F 
which characterizes the inverse kinetic equation results uniquely determined. Furthermore, 
we have proven the positivity of the kinetic distribution function. An open issue remains, 
however, whether the inverse kinetic equation ([2]) satisfies an H-theorem, i.e., the entropy 
results a mononically increasing function of time. Finally, as an additional development, we 
have shown that, consistently with the assumption that the fluid fields are strong solutions of 
INSE, the mean-field force can be expressed in such a way to satisfy explicitly also the energy 
equation. The result appears significant from the mathematical viewpoint, the physical 
interpretation of the theory and potential applications to the investigation of complex fluids, 
such as for example those treated in g], 10, llll)- In fact, it proves that the inverse kinetic 
theory developed in pj] can be given an unique form which applies to classical solutions of 
INSE. 
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